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In this paper, we deal with the problem of implementing an abstract machine for a stochastic ver- 
sion of the Brane Calculus. Instead of defining an ad hoc abstract machine, we consider the generic 
stochastic abstract machine introduced by Lakin, Pauleve and Phillips. The nested structure of mem- 
branes is flattened into a set of species where the hierarchical structure is represented by means of 
names. In order to reduce the overhead introduced by this encoding, we modify the machine by 
adding a copy-on-write optimization strategy. We prove that this implementation is adequate with re- 
spect to the stochastic structural operational semantics recently given for the Brane Calculus. These 
techniques can be ported also to other stochastic calculi dealing with nested structures. 

1 Introduction 

A fundamental issue in Systems Biology is modelling the membrane interaction machinery. Several 
models have been proposed in the literature ifTTl [141 13[ ; among them, the Brane Calculus [4] has been 
arisen as a good model focusing on abstract membrane interactions, still being sound with respect to 
biological constraints (e.g. bitonality). In this calculus, a process represents a system of nested com- 
partments, where active components are on membranes, not inside them. This reflects the biological 
evidence that functional molecules (proteins) are embedded in membranes, with consistent orientation. 

In the original definition of the Brane Calculus |4] (which we will recall in Section[2]) membranes in- 
teract according to three basic reaction rules corresponding to phagocytosis, endo/exocytosis, and pinocy- 
tosis. However, this semantics does not take into account quantitative aspects, like stochastic distribu- 
tions, which are important for, e.g., implementing stochastic simulations. 

A stochastic semantics for the Brane Calculus has been provided in [2], following an approach pio- 
neered in lH (but see also [8, 12] for Markov processes). Instead of giving a stochastic version P Q 
of the reaction relation, in this semantics each process is given a measure of the stochastic distribution 
of the possible outcomes. More precisely, we define a relation P — > pt associating to a process P an 
action-indexed family of measures pt: for an action a, the measure \i a specifies for each measurable set 
S of processes, the rate /J. a {S) G M + of a-transitions from P to (elements of) S. An advantage of this 
approach is that we can apply results from measure theory for solving otherwise difficult issues, like 
instance-counting problems; moreover, process measures are defined compositionally, and in fact the 
relation P — > /J. can be characterized by means of a set of rules in a GSOS-like format. We will recall this 
stochastic semantics and its main properties in Section[3] 

In this paper, we use this new semantics for defining a stochastic abstract machine for the Brane 
Calculus, so that it can be effectively used for in silico simulations of membrane systems. Defining an 
ad hoc abstract machine for the Brane Calculus would be a complex task; instead, we take advantage of 
the generic abstract machine for stochastic process calculi (GSAM for short) introduced in 1 131 [TOl as 
a general tool for simulating a broad range of calculi. This machine can be instantiated to a particular 
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calculus by defining a function for transforming a process of the calculus to a set of species, and another 
for computing the set of possible reactions between species. 

An important aspect is that this abstract machine does not have a native notion of compartment, 
which is central in the Brane Calculus (as in any other model of membranes). To overcome this problem, 
we adopt a "flat" representation of membrane systems, used also in iflOl . where the hierarchical structure 
is represented by means of names: each name represents a compartment, and each species is labelled 
with the name of the compartment where it is located, and the name of its inner compartment (if any). 
So names and species are the nodes and the arcs of the tree, respectively. This technique can be used for 
representing any system with a tree-like structure of compartments. 

However, this approach does not scale well, as the population of species may grow enormously: for 
instance, a population of n identical cells would lead to n species, all differing only for the name of its 
inner compartment, instead of a single specie with multiplicity n. For circumventing this problem, in 
Section [4] we introduce a variant of the GSAM with a copy-on-write optimization strategy — hence called 
COWGSAM. The idea is to keep a single copy of each species, with its multiplicity; when a reaction has 
to be applied, fresh copies of the compartments involved are generated on-the-fly, and reactions and rates 
are updated accordingly. In this way, the hierarchical structure is unfolded only if and when needed. 

In Section [5] we show how the Brane Calculus can be represented in the COWGSAM, and we will 
prove that the abstract machine obtained in this way is adequate with respect to the stochastic semantics 
of the Brane Calculus; in this proof, we take advantage of the compositional definition of this semantics. 

Conclusions and final remarks are in Section [6] 

2 Brane Calculus 

In this section we recall Cardelli's Brane Calculus [4] focusing on its basic version (without communi- 
cation primitives, complexes and replication). 

First, let us fix the notation we will use hereafter. Let S be a set of sorts (or "types"), ranged over 
by s,t, and T a set of 5-sorted terms; for t G S, T t C T denotes the set of terms of sort t. For A a set of 
symbols, A* denotes the set of finite words (or lists) over A, and (ai,...,a n ) denotes a word in A*. For a 
word (h,.. .,t n ) in S*, we define T(t u .. 4n ) =T h x---x T, n . 

Syntax The sorts and the set B of terms of Brane Calculus are the following: 

t ::= sys | mem 
ct,t::=0 | g\x | »„.a | »^(r).a | s„.a | Q^.a | ® n (-r).a 
P,Q::=o | P°Q | adPD 

The subscripted names n are taken from a countable set A. By convention we shall use M,N, ... to 
denote generic Brane Calculus terms in B. 

A membrane can be either the empty membrane 0, or the parallel composition of two membranes 
a\x, or the action-prefixed membrane e.G. Actions are: phagocytosis », exocytosis Q, and pinocytosis 
©. Each action but pinocytosis comes with a matching co-action, indicated by the superscript J -. 

A system can be either the empty system o, or the parallel composition P Q, or the system nested 
within a membrane a^Pl). Notice that, differently from [4], pino actions are indexed by names in A. 
In H , names are meant only to pair up an action with its corresponding co-action, hence a pino action 
does not need to be indexed by any name. Actually, names can be thought of as an abstract representation 



Sorts ::S 
Membranes : : B 
Systems :: B sys 
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; — : : : - : ; : ; — — (red-pha£o) 

»i(p).T|T dgl) v„.o\o <lPD — ► r\r dp<la\ao<lPDD ° 2D 

(rcd~cxo^ 

»^.T|T <l8 B .(T|ob<lPD ° 2D — ► a|a |T|T (l2D °P 

(red-pino) — - — - — — (red-loc) 



@(p).a|a <1/ 5 D^a|ao(!p(!*DoPD ' adPD^adQD 
P^Q P = P' P'^Q' Q' = Q 

(red-comp) — — (red-equiv) 



PoR^QoR P^Q 



Table 1: Reduction semantics for the Brane Calculus. 



of particular protein conformational shapes; hence, each name can correspond to a different biological 
behaviour. Therefore, if we want to observe also kinetic properties of processes, it is important to keep 
track of names in pino actions. 

Terms can be rearranged according to a structural congruence relation; the intended meaning is 
that two congruent terms actually denote the same system. Structural congruence = is the smallest 
equivalence relation over B which satisfies the axioms and rules listed below. 

P°Q = Q°P P° (Q°R) = (PoQ)or poo=p 
o\t = t\o c|(t|p) = (cj|t)|p a|0 = cr 
«.„ k P = Q o = t P = Q a = t 

oaoD = o — — j — — - — — — 

Por = qor <j|p = T|p a«PD = T(l2D 

a g {s, 7 ,s n ,s,t}„ £ A a = t p g {»^,©„}„ GA p = v a = T 
a.o = a.z j8(p).cj = j8(v).t 

Differently from @, we allow to rearrange also the sub-membranes contained in co-phago and pino 
actions (by means of the last inference rule above). 



Reduction Semantics The dynamic behaviour of Brane Calculus is specified by means of a reduction 
semantics, defined over a reduction relation ("reaction") ^ C B sys x B sys , whose rules are listed in 
Table [T] Notice that the presence of (red-phago/exo/pino) and (red-equiv) makes this not a structural 
presentation, since these rules are not primitive recursive in the syntax (i.e., structural recursive) as 
required by the SOS format. 



3 Stochastic Structural Operational Semantics for the Brane Calculus 

In this section we recall the stochastic structural operational semantics for the Brane Calculus, as defined 
in 0- Following lH, we replace the classic "pointwise" rules of the form P P 1 with rules of the form 
P — > /J,, where p is an indexed class of measures on the measurable space of processes. We assume the 
reader to be familiar with basic notions from measure theory; for a brief summary, see Appendix [A] 
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a G {»„,»„, K>£ | n G A} 

(zero) : — : (pref) 



''mem ®mem QJ.CT > me m [Of] 

j3 6 {»^,©„ I n £ A} . g -^mem /*' T ^ mem \l" 

j8(T).a ^mem [j8]J ff| T -> me m /i'cr©T JU" 

C ^mem V P 7- S y S 

^sys M 



( vold ) ( loc ) — „. ^ r~7 ^ » — (comp) 



^sys 



Wsys aCPD ^ sysi U@pV P°Q->sysV'p®QH" 



Table 2: Stochastic structural operational semantics for Brane Calculus 



The set of action labels for the Brane Calculus will be denoted by A and can be partitioned with 
respect to the source sort (i.e., either systems or membranes), as follows: 

^sys — {id : sys — > (sys)} U {ph n : sys — > (sys, sys) | n G A} U 
{ph,f : sys — > (mem, mem, sys, sys) | n G A} U 
{ex„ : sys — > (mem, sys, sys) | n G A} 

Amem - {»«, »«, : mem ~> (mem) | n G A} U 
{s^", @„ : mem — >■ (mem, mem) | n G A} 

Let a range over A, and arid) denote its arity. To ease the reading in the following we will use the 
notation A a (T,L) to denote the set of measures A(7( fl ^,®" =1 E ?1 ), for ar{a) = t — > (t\, . . . ,t n ). 

Let B/= be the set of =-equivalence classes on B. For M G B, we denote by [M]= the =-equivalence 
class of M (sometimes dropping the equivalence symbol when clear from the context). 

Definition 3.1 (Measurable space of terms). The measurable space of terms (B,IT) is given by the 
measurable space over B where II is the O -algebra generated by B/=. 

Notice that B/= is a denumerable partition of B, hence it is a base (a generator such that all its 
elements are disjoint) for IT. Any element of IT can be obtained by a countable union of elements of 
the base, i.e., for all M G IT there exist {M,} lG /, for some countable /, such that M = |J ;g/ [M,]=. As a 
consequence, in order to generate the whole IT we can simply compute all these unions, without the need 
of any closure by complement. 

A similar argument holds for the product space (B^ f „\,(8>" = in f; ), where t[ G {mem, sys} (1 < 
i < n)\ indeed (g>" =1 n, ; can be generated from the base M^^/^ ^, where = (fl jfn) C B<, lr .. A ) x 
%,...A> is denned by 

(M b . . . ,M n ) = {h ,..., tn) (N u ... ,N n ) iff Mi = N h for all 1 < i < n , 

which can be easily checked to be an equivalence relation. =/ h ^-equivalence classes are rectangles, 
i.e. [(Mi , . . . ,M„)]= ( tn) = [M\]= x • • • x [M n ]=, therefore the product measure (g)" =1 Tl ti is well defined. 
For sake of simplicity in the following we write [(Mi, . . . ,M„)]= in place of [(Mi, . . . ,M n )]= { tn) , and 

% ; „>/- in P laceofB ^,..A)/- (fl ,^„)- 

The operational semantics associates with each membrane a family of measures in A Amem (B,n), 
and with each system a family of measures in A As y s (B,IT). This can be represented by two relations 
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-^ mem : T mem -> A A — (B,n), -^ sys : r sys ->■ A A ^(B,n), defined by the SOS rules listed in Table|2] (In 
the following, for sake of readability, we will drop the indexes mem ,sys)- In these rules we use some 
constants and operations over indexed families of measures, that we define next. For a set (of labels) A, 
let us denote by A A (B,n) the set FLeA A a (B,n) of A-indexed families of measures over (B,n). Given a 
family of measures JJ. G A A (B,IT) and a G A, the a-component of /J. will be denoted as \i a G A a (B,n). 

Null: Let (O mem G A Amem (B,n) be the constantly zero measure, i.e., for all aGA mem such that ar(a) = 
t-+{ti,...,tn) and M e <S)"=i n tr (co mem ) a (M)=0. 

Prefix: For arbitrary n G A, a G {s,^,^- 1 }, and j3 G ©}, let the constants [a„] a , \fin\o G A Amem (B,n) 
be defined, for arbitrary X,Y G B mem /=, as follows: 



Parallel: For jU,ju' G A Amem (B,n), let jU ct (D t ju' G A Amem (B,n) be defined, for n G A, a G {s,^,^- 1 }, 
j8 G {a- 1 ,®}, andX,y GB mem /^, as follows (where for allX,T: Xi T = U{M= I ct|tgX}): 



Void: Let 0) sys G A As y s (B,n) be defined by ((O sys ) a (!M) = for any a G A sys , such that ar(a) = t — > 
,r n ) f and fWG(g)? =1 n fj . 

Nesting: For v G A A ™ m (B,n) and jU G A A ^(B,n), let n@%v G A A ^(B,n) be defined, for X,Y G 
B mem /= and Z, W G B sys /=, as follows: 




(la n } a ) Pm (XxY) 












(M ff ® T ju') «„ (X) = AV, (X, T ) + ^ (X, a ) 




Ui@?v)id(X) = Hainan) + 



v @ „(X'xX") + 



X'3A:"(![«] s Do[P] s D=X 




i(n) 



x'|A:"ay"Doy'=z 
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Composition: For £ A A ^(B,n), let \lp®Q \l' G A A ^(B,n) be defined, for X,Y G B mem /= and 
Z,W e B sys /^, as follows (where for all W, Q, W oQ = \J{[P}= \P°Q£ W}): 

(Atp0 e At') P h„ (ZxW) = Ai P h„ (Z x Wo 6 ) + Mph„ (Z x Wop) 

(jUp^q n') pb x (XxY xZxW) = Ai ph ± (X x Y x Z x W oQ ) + ^ (XxYxZx Wo P ) 

(Atp^Q At')ex„ (X X Z X W) = Aiex„ (X X Z X Woq) + jU^ (X X Z X Wop) 

(H P ® Q H') id (X) = Hid(Xo Q ) +n' id (Xo P ) + 

ne\ n pK (Yi x Y 2 ) • tt ' (Xi x X 2 x Zi x Z 2 ) 

— h 

l(n) 

n ^> Mph^ (Xi x x 2 x y t x y 2 ) • ^ K (z x x z 2 ) 

^ l(n) 

These operators have nice algebraic properties (e.g., ii' ff (D T \l" = \l" ' T (Dcr A 4 '* (/j'ctCDt M")a|T®p A 4 "' = 
M'ct^tIp (A 4 '' VDp A 4 "')' • • • )» an d respect the structural congruence (e.g., if P = P f and Q = Q f then jl'p®Q 
\l" = il'pi®@ /J."). We refer to [2] for further details about these properties, very useful in calculations. 

The next lemmata state that the stochastic transition relation — > (and hence the operational semantics) 
is well-defined and consistent, that is, for each process we have exactly one family of measures of its 
continuations, and this family respects structural congruence. 

Lemma 3.2 (Uniqueness). For a G A such that aria) = t — > (ti, . . . , t n ), andM G M t , there exists a unique 
At G A A ' (B, n) such that M — >■ /x. 

Lemma 3.3. IfM = N and M ->■ fx, then N ->■ A 4 - 

This operational semantics can be used to define the " traditional" pointwise semantics: 

M^(M u ...,M n ) <^M^pl and Al a ([(M 1 ,...,M„)] s ) = r 

and it is conservative with respect to the non-stochastic reduction semantics. 
Proposition 3.4. For all systems P,QG B mem , ifP^ji and \i id ( [Q] ) > ?/z<?« P^ Q- 



4 The COW Generic Stochastic Abstract Machine 

In this section we present a variant of the generic stochastic abstract machine (GSAM), oriented to 
systems with nested compartments. 

The GSAM has been introduced in lfT3l[T0l for simulating a broad range of process calculi with an ar- 
bitrary reaction-based simulation algorithm. Although it does not have a native notion of "compartment", 
nested systems can represented by "flattening" all compartments and species into a single multiset, where 
each species is tagged with names representing their position in the hierarchy, as shown in ifTUl . The idea 
is to represent each compartment as a different species, keeping track of their position in the hierarchy 
by means of (node) names. These names are ranged over by x,y,z, ■■ • , and are different from names in 
actions. As an example, a system adTdDD is represented as the multiset {(7<0 h-> 1 , Trd^D h-> 1}, which 
means "there is one cell with membrane a located in the compartment x and whose compartment is y, 
and one cell with membrane z positioned in y and whose compartment is z". Reactions can happen only 
if the names tagging the involved species match according to the required nesting structure. 
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Unfortunately, this approach does not scale well, as the population of species grows. Let us consider 
a system composed of n copies of the same cell, e.g., n • (crflD) (where n can be easily in the order of 
10 3 -10 4 ). In the original GSAM idea, this should be represented in the machine as a single species with 
multiplicity n, and each possible reaction is represented once but with propensity given by the law of mass 
action taking into account the species' multiplicity n. Instead, the "flat encoding" above yields n different 
species cfl* D, . . . , (7<K D, each with multiplicity 1; the set of reactions explodes correspondingly. 

For circumventing this problem we introduce a variant of the GSAM with a copy-on-write strategy — 
hence it is called COWGSAM. The idea is to keep a single copy of each species, with its multiplicity; the 
same applies to reactions. When a reaction has to be applied, the compartments involved are "unfolded", 
i.e., fresh copies of the compartments are generated and the reaction set is modified accordingly; then, 
the reaction can be applied. In this way, the hierarchical structure is unfolded only if and when needed. 

In order to implement this idea, we have to modify the notion of machine term, reaction and reaction 
rule. The COWGSAM (with the Next Reaction method) is shown in Figure [T] 

First, for generating fresh names, we have to keep track of those already allocated. To this end we 
introduce environments, which are finite sets of names. Then, the machine state is represented by a 
machine term T, i.e. a quadruple E h (t,S,R) where E is an environment; t is the current time; S is a 
finite function mapping each species / to its population S(I) (if not null); and R maps each reaction O 
to its activity A, which is used to compute the next reaction. (Notice that the syntax of species / is left 
unspecified, as it depends on the specific process calculus one has to implement.) We say that a machine 
term E\- (t,S,R) is well-formed if for all x t ,Xj G E :xt = xj i = j, and the free names in S,R appear in 
E. In the following, we assume that machine terms are well formed. 

Each reaction is a quadruple (Si , r,f, S2), basically representing a reaction Si — > S2, where 

• Si and S2 denote the reactant population and product population respectively; 

• r denotes the rate (in s~ l ) of the reaction; 

• / is a function mapping machine terms to machine terms; this functions implements any global 
update of the machine term after the reaction (if needed). 

The transitions of the abstract machine are represented by a relation T T between machine 
terms, indexed by the propensity a and reaction rules. This should be read as "T goes to T' with rate 
a, using the rule O". This relation is defined by (Reaction rulei in Figure [I] where the function next(T) 



selects the next reaction, i.e. it returns a pair (0,t') where O = (S\,r,f,S2) is the reaction to happen 
first among all possible reactions in T, and t' is the new time of the system. Once the reaction has 
been selected, we have first to create the separate (private) copies of the compartments involved, and to 
update the reaction set accordingly. This is done by the functions cow(_) and dup{ _), which implement 
a deep copy-on-write: cow(E h (t,S,R),Si) is a machine term E' h (t,S',R') representing the same state 
as E h (t,S,R), but in 5' the species indicated in Si are unfolded; E' contains all names which have 
been generated in the process, and R' is the new set of reactions. (Actually, the freshly generated copies 
represent the instances which are not involved by the reaction; this simplifies the reaction application.) 
An example of the action of cow(S) is depicted in Figure [2] 

At this point, the reaction is executed, by removing the reactants Si from the machine term (via 
the operation Q), adding the products S2 (via the operation ©) and updating the current time of the 
machine. The function / performs any global "clean-up" and restructuring to the machine term that 
may be required by the reaction (e.g., garbage collection, elimination of names not used anymore,. . .). 
Moreover, since a reaction can rearrange the hierarchy structure, possibly creating new compartments 
and deleting others, we have to add to the environment any fresh name introduced in the products. 
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T 
E 
S 
R 
O 



=E\-(t,S,R) 
-X\ , . . . ,x n 

={h ^h,...,I N ^i N } 
={0i \->Ai,...,0 N i->A N } 
=(Si,r,f,S 2 ) 



((Si,r,f,S 2 ),a,t')=next(t,S,R) E' h (t,S',R') = cow(E h 



(Machine term) 
(Environments) 
(Populations) 
(Reactions) 
(Reaction) 

(Reaction rule) 



cow 



next(t,S,R) = {0,a,t') if R(0) = (t',a) and t' = mm{t \ R{0) = (t,a)} 

cow(E h (t,S,R),Q) = Eh (t,S,R) 
{Eh {t,S,R),{pd y z D^ j}USi) = cow(E'ufn(S") h (f,S" Uinit(L' UL" ,(t,S' ',R)),Si) where 

if 5((7dyD) > 1 for some a,x : 
let i = S((7d*D) and/ (£E 

(S',R r ) = dup{E h (f,S{crd*I) ^ l,crd*,D ^ i- \},R),y,y') 

L' = react ions (o dyD i-> /— 1,5) 
otherwise let (5',/?') = (S,R),L' = in 
letn = S(pdp),£' = ,EU/n(S'), andz g £' 

(5",/?") = dup(E' h (?,5'{p(ip m- j,p(P,D m- n-y},/?'),^) 

L" = reactions(pd y z ,b \->n — j,S') 

dup(Eh {t,S,R),y,y') = {S' U S" ,RU R" Uinit(L,(t,S' ,R))) 

where S' = SU {pd^D ^ i | S(ad^) > 0,S(pd£,D) = ?},w' £ E, 

L = react ions (pd^,D i-> z',S ), 

(5", /?") = rfiip (£ U {w'} h (f , 5, /?) , w, w') 

Figure 1 : The COW Generic Stochastic Abstract Machine, with the Next Reaction method. 

Finally, the term can be "normalized" by collapsing equivalent copies of the same subtree into a 
single copy (with the right multiplicity), by the function norm{_). In its simplest form, norm(_) can be 
the identity, i.e., no normalization is performed at all. Although this is correct, it can lead to unnecessary 
copies of the same subtrees. We can define norm(_) such that 

norm({h i-> i 1: I 2 ^ 12} US,/?) = norm({h ^ h + i 2 }U S,R[I 2 /h}) if h = h 

where the equivalence between species can be implemented by comparing the subtrees starting from I\ , 
I 2 , e.g., by calculating a suitable hash value. We leave this refinement as future development. 

In order to implement the Next Reaction method, each reaction O is associated with a pair R(O) = 
(a,t), where a is the reaction propensity and t is the time at which the reaction is supposed to occur. The 
function delay(r, a) computes a time interval from a random variable with rate r and propensity a. 
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{h !->• h,..., i N ^i N } ®(t,s,R) =h^ke...ei N ^i N e(t,s,R) 



I^i®(t,S,R) = < 



' (t,S',RUupdates(I, {t,S',R))) if S(I) = i and S' = S{I ^ i' + i} 
(t,S',RUinit(L, (t,S',R))) if / g dom(S),S' = S{I ^ i} 

and L = reactions(I i-> z',S) 



0, 5,7?) {/i i-> /i,..., I N ^i N } = (t,S,R)Qh^hQ...Ql N ^i N 
(t,S,R)Ql^i= (t,S',R(Jupdates(I, (t,S',R))) 

if 5(7) = i', i' > i and S' = S{I ^ /' - i} 

init(L, (t,S,R)) = {04 | O G L and a = propensity(0,S) and 

= (S\,r,f,S 2 ) and f' = t + delay(r,a)} 

updates(I, (t,S,R)) = {O ^ (t',a) \ R(0) = (t",a) and O = (S u r,f,S 2 ) and 5i(7) > and 

if t" > t then a = pwpensity(0,S) and t' = t + (a/a')(t" - 1) 
if t" = t then a 1 = propensity(0,S) and t' = t + delay(r,a)} 



pmpensity((Si,r,f,S 2 ),S) = r- 



S*(h 
ji 



S*(I n ) 
jn 



if Si ={h \->ji,...,I n ^ jn} 



s -W) = |s(0<W 



if x = root 



y 

S(Gdp)-S*(pd z x D) if jc i-root and 5(pdp) >0 



Figure 1: The COW Generic Stochastic Abstract Machine (cont.). 




Figure 2: Example copy-on- write of a subterm. 



This general structure can be instantiated with a given process calculus, just by providing the defini- 
tion for the missing parts. Given a set Proc of processes, we have to define: 

1. the syntax of the species I (which may be different from that of processes); 

2. a function speciesE, x {P) mapping a process P 6 Proc to a species set located in x; 

3. a function reactionsil >->■ i,S) for computing the multiset of reactions between a (new) species I 
with multiplicity i and a population of (existing) species S. 
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The function species is used to initialise the abstract machine at the beginning of a simulation. If we aim 
to simulate the execution of a process P G Proc, the corresponding initial state (rooted in root) is 

<\P\)root = fn(J) h J® (0,0,0) where / = species®. root (P). 

The reactions function is used to adjust the set of possible reactions dynamically. 

5 Implementing the Stochastic Brane Calculus in COWGSAM 

In this section, we show how the COW Generic Stochastic Abstract Machine can be used to implement 
the Stochastic Brane Calculus, following the protocol described in the previous section. 

5.1 Encoding of the Stochastic Brane Calculus 

Syntax of species We define the species for the brane calculus, which in turn lead us to introduce 
complexes and actions. Notice that (despite the deceiving syntax) species are not systems and complexes 
are not membranes; nevertheless, actions are a subset of membranes. 

I:- -Cd x y l> (Species) 
C ::= A\, . . . ,A n (Complexes) 
A ::= s„.(7 | »^(r).a I Q„.a | »^.a | ©„(r).a (Actions) 

Node names can appear in the species in S and in reactions R. 

The species function We can now provide the definition of the translation of a process P G B sys into a 
species set. Basically, each compartment is assigned a different, fresh node name; therefore, the function 
speciesE, x (P) is parametric in the set E of allocated node names, and the name x G E to be used as the 
location of the system P. The name x changes as we descend the compartment hierarchy. 

In order to capture correctly the multiplicity of each species, we assume that systems are in normal 
form. Basically, this form is a shorthand for products where n copies of the same system, i.e., Po-oP, 
are represented as n ■ P. 

Normal Systems :: B^ ys Q ::= m ■ o^QiD ° . . . ° n k ■ a k dQ k l> 

where k>0 and for i^j: oAQfi ^ OyflQyD 

A system in normal form can be translated into a system just by unfolding the products. For Q G Bs ys a 
system in normal form, let \Q] G B sys defined as follows: 

n times 

[o] = o \ n - ade'D o Q"] = a« \Q'] D o • • • o ad \Q'] D o \Q"] 

Proposition 5.1. For all P G B sys , there exists a system in normal form Q G B" ys such that \Q] = P. 
As a consequence, we can give the definition of species on systems in normal form, as follows: 

speciesE } x( ) — 

s pedes Ejc(n- odQiboQ 2 ) = {s(a)fl*D !->■ n} U species EW { y y iy (Qi) U speciesE^{Qi) 

with y g E and fn(species E ^{y} : y{Q\)) nfn(species E , x (Q2)) C {x} 
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The condition in the second case ensures that two different compartments are never given the same 
name — any name clash has to be resolved by an a-conversion. The function converts a membrane 
into a set of complexes: 

*(O) = s(o\o') = s(o)Us(o') s{n.o) = {n.o} 

The reactions function The next step is to define the function reactions(I >->■ i,S), for I a species with 
multiplicity i and S a population. 

reactions s(h i,S) — unary e(I\) U binary e(I\ , S) 

unary E (h) 4 {({I, ^ \},r n ,id,{(s(o)UU[)dp ^ l,^(p)^D ^ 1}) | 

h=U^p, U X = {®„(p).(j}UU[,w^E} 
binary E (h,S) 4 {({/j ^ l,/ 2 ^ l},r„,/, {(s(z) U I/( Uj((t) U U 2 )« x y D ^ 1}) | 7 2 G <fom(S), 

(/! = c/idjD,^ = c/ 2 ^D) v (/ 2 = = u 2 «yj), 

Ui = {^.%}UU[,U 2 = {^ n .(j}UU 2 ,f = XT.T[w:=x]} W 
{({/i i-> l,/ 2 ^ l},r„,i'rf, 

{(^Uf/O^ -> MpW -> l,(^(a) U^)ap M- 1}) | 7 2 G <fom(S), 

(/j = c/ia*D,/ 2 = c/ 2 dp) v (i 2 = cMJMi = u 2 i x z t>), 

C/ 1 = {^(p).T}Uf/ 1 ',f/ 2 = {»„.<7}Uf/^ w(£E} 

In the case of Brane Calculus, the unary reactions are only those arising from pinocytosis, while binary 
reactions arise from exocytosis and phagocytosis. In both cases, the multiplicity of each reactant is 1, so 
the multiplicity of 1\ is not relevant. Exocytosis merges two compartments; this is reflected by the fact 
that the "rearranging" function / substitutes every occurrence of the name w in T with x. On the other 
hand, pinocytosis and phagocytosis create new compartments; to represent the new structure, we choose 
a fresh name w representing the new intermediate nesting level, and reconnect the various compartments 
accordingly. Therefore, for any reaction (Si,r,f,S 2 ) G reactionssih | -> i,S), fn(S 2 ) \E is either (in 
the case of exocytosis) or the singleton {w}. 

5.2 Adequacy results 

Before proving the correctness of our implementation, we have to define how to translate a species set 
back to a system of the brane calculus. 

Let S be a non empty species set. A root name of S, denoted by root(S), is a name x such that 
S(Cd*l)) > for some C,y, and for all z,C : C'4 Z X D g dom(S). The next result states that root(.) is well 
defined on the species sets we encounter during a simulation. 

Lemma 5.2. For all P G B sys : 

1. ifP^»: root(<\P\) x )=x. 

2. if(\P\)x E h (t,S,R) andS^ 0, then root(S) = x. 

Proof. (1.) is trivial by definition. (2.) It is enough to check that the reaction rules introduced by 
reactions^ do not change the name of the root, nor introduce new ones. □ 
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We can now define a function [_] which maps complexes to membranes, and a function mapping 
species sets to systems; the latter is parametric in the name x of the root of the system. 

lA h ... ,A„] =Ai| ... \A n {Sj x 4 J] 5(C(ip) • ([CM%D) 

CCpG(/om(S) 

where the notation n ■ P is a shorthand for P . . . P, n times. 
Lemma 5.3. For all P € B sys : 

1. UpUx=p- 

2. if<\P\) x ^Eh (t,S,R) then {SJ X is well defined. 

Proof. (1.) is easy. (2.) It is enough to check that the reaction rules introduced by reactions{_) do not 
introduce loops (i.e., the order among names is well founded). □ 

We can now state and prove the main results of this section. 

Proposition 5.4 (Soundness). For all P G B sys , §P\ X — > E h (t,S,R) then there exists n such that 
P^H and Hu([[Sl x ]) = a. 

Proof. The proof is by cases on which reaction rule O is selected by the function next. By additivity of 
measures, we can restrict ourselves to when the whole process P is the redex of the reduction. Let us see 
here the case when P is the redex of a (red-pino) (another is in Appendix [B]). 

Let P = © n (p).o\(Jo(lP'V and let us assume that Ob does not exhibit a ®„ action. Then, the translation 

of P is t\P\) x = E h species 9jX {P) © (0, 0, 0), where 

species® iX (P) = {s(@ n (p).o\o )<l x y V h-» 1}U species ^^{P 1 ) 

= {(s(@ n (p).o)Us(o Q ))(ip ' ^ \}U species {yU {P') 
= {({® n (p).a}Us(a ))(Lp h-> 1} Us pedes {yhx (P r ) 

Let h = ({®„(p).o}l)s(o ))<Lp and J P = species {y}x (P'); then <\P\) X =h^l®J P r® (0,0,0) = J P © 
(0,S',R') where 

5' = {/i i ^ 1} 

R' = init(L,(Q,S' ,%)) = {0 L ^ (t x ,a{)} 

L = ({h ^ l},r n ,id,{(s(o)Us(o ))dp ^ l,s(p)« y J k> 1}) 
L = reactions(h ^ 1,0) = unary(l\) = {{{h ^ l},r„,irf,{(j((T) U j(ob))<IyD i-> 1, j(p)d^D h-> 1})} 

Now, the reaction in ^ * s ®l (otherwise it would involve P', not the pino of the whole P). 



This means that <\P\) X T is derived by means of an application of the (|Reaction rule)) as follows, where 



Si = {h ^ 1} andS 2 = {(j(ct)U *(ob))<ip ^ M(P)0 ^ 1}: 

( (5 b r H ,id,S 2 ), a h h) =next(0,S',R') (E' h (0,S",#")) = cow(£ h (0,S',/?'),5i) 

g r-{-Me(o,y,/?') ai '^' r " aA )) wom (£u/n(5 2 )h ({/ P /}®s 2 e((fi,^,^)e^i))) 

where 5" = 5' and /?" = /?'. 



94 



Implementing Stochastics Brane Calculus in COWGSAM 



{Jp,}®S 2 (B({tuS"X)QS l ) = 

= Jp,(B{(s(o)Us(o Q ))(ip^l,s(p)(liD ^ l}(B((t u S',R')Q{h 1}) 

= J P © {{s{o) U j(ob))<IJD ^ M(p)0 >-> 1} (/i , {/i ^ 0}, {<9 L i— > (t 2 ,a 2 )}) 

= Jp'®J'®(tu{h^0}:{O L ^(t2,a 2 )}) 

Now let us define Q as 2 = alobtfpfloD QpD, then (|0^ = species^ x (Q) © (0,0,0) where 

species<D )X (Q) = species<6, x (c\CQ<ip(l°]) °P'D) 

= {s(ff|a )<IyD h-> l}UjpedeJ Wi y(p(I«D °P') 

= {(s(a)Us(cb))fl*D H> l}UjpedeJ W)> (p(I*D)Uj/»edeJ Wj> ,(F') 

= {(5(a)U5(a ))dyD !-»■ l}U{j(p)(ip i-> l}U0U/ P /=/u/ F 

and hence Q = {J P © / © (t\ , {l\ H> 0},{C?l h-> (t2,aa)Y)\ x . It remains to prove that r„ = Hm([Q\). Let 
us notice that the derivation of P — > p is actually as follows: 

o„(p).(T->- [©„]§ a ^p" 

(par) 



©„(p).a|a -> [@/a e „(p). 5 (Da„M" P'^ju' 

(l° c ) — i — -r~ 

®„(p).a|a (!P , D->p 

W herep = p'@p [p) ala \[@XMP)^o o n- Then: 

^([alaodpdoDop'D]) = (p'@®^ ff l ff0 ([@X @n(p ). CT ©aoM"))^(^l^o(Ip^D°P'D]) 

= ^([^kodp^D°P / D]) + ([@n]& @n (p). ff ®cT M") @n (^ko] x [p]) 
= ([•»]So,(p).ff®a b M")o.([ff|ob] x [p]) 

= ([® n ] P o)® n ([o\Go} >< [P])+<,(^l^o] X [p]) 

= r n + ^ n (i a \ a o\ x = r « 

where the last equivalence holds because ju@ ([cr|<7o] x [p]) = because we assumed that the reaction 
does not involve Oq. □ 

Proposition 5.5 (Progress). For all processes P,Q, ifP^Q then there exists a reaction O and a term 
T such that (\P\) X ^ T and Q = {Tj x . 

Proof. By induction on the derivation of P ^ Q. Let us see the case of (red-pin), the others being similar. 
LetP = ®(p).o\o dP'\) and Q = o\o <Lp<loD°P'D. Then, 

<\P\, X = species {yU {P') © (0,{({o„(p).ff} Us(a ))<ip ^ 1},{0 L ^ 

whe re Q L = ({({® „(p).a} Us(a ))(l*i) ^ l},r„,id,{(s(o)Us(oo))<l*D ^ M(p)Qh> 1}). Then, by 
the ((Reaction ruleb we can take J = species^ (P') © (fi, {0(a) Us(a ))(!;*D h-> l,j(p)CD h-> 1},0). It 
is easy to check that Q = lT} x . □ 

Proposition 5.6 (Completeness). For all processes P,Q, if P — > p and p,y([<2]) > then for all node 
name x, there exists a reaction O and a term T such that t\P\) x T, Q= fTj x and a = Pid{[Q})- 



Proof If P -> p and p id ([Q]) > then P — ► Q by Prop. [^4] By Prop. ||3[ we have that for some a,0,T, 
(|PD* ^> r and Q = [r]*. But then a = p,y ([£>]) by soundness (Prop. 



5.4 1. □ 



M. Miculan and I. Sambarino 



95 



5.3 Example 

We conclude this section with an example. Let P = 10000- »„.S m d» fe (IDD ° 100 - ((»^(Q^) | S J -)a» i: (IDD). 
Then, its reductions in the Brane Calculus are as follows: 

10000 ■ v n .v m 4v k <Ln o 100 -{(K{Vm) I s^MDD) 

—►9999 • s„.Q m ds fc <IDD °99 • ((»^(s^) | s- L )<3» jt <IDD) sM^flSmfl^dDDD s^dDD 
—►9999 • s„.Q m (Is fc <IDD °99 • ((s^(s^) | o^ds^ O 555 ) ° Q " L(I2 ■ °d*I>D 

The translation of P is = E h species®, x (P) (0,0,0), where 

species^P) = species {x} (10000 ■ S^.^ds^DD o 100- (» x (s,t) | O^aSfcdDD) 

= {*(»„.a m )<I*D ' ^ 10000} u^eci« {lJ(}o ,(Sfc(ID)u{j(si(o^) | s x )<I*I> ^ 100} 

U spec iesr XiZ \ t z ( C D ) 
= {{»„.s m }(I^ i ^ 10000, {« x (s x ), o x }fl*I> m- 100, m- 1, 

{»,}a z v D^i} 

Let /j = {» n .S m }C*D, I 2 = O x }«P, / 3 = {®*}0, /4 = {»*}« Z V I>, r„ = Ids" 1 , = 5s' 1 and 

r m = 5s~ 1 ; then 

(\p\) x = h ^ iooooe/ 2 ^ iooe/ 3 ^ ie/ 4 ^ 10(0,0,0) 
= h^>- iooe/ 3 1-> i©/ 4 i-> i©(o,Si,/?i) 
= / 3 h+ i©/ 4 ^ ie(o,5 2 ,/? 2 ) 

= / 4 ^le(0,5 3 ,/?3) = (0,54,/?4) 

where 

Li = reactions(h >-)• 10000,0) = 

51 = 5{/i ^ 10000} 

/?i =/JUi'/iif(Li,(O,5i,0)) = 

L 2 = reactions^ i-> 100, Si) 

= ({/ 2 m- l,/i h- 1}, 10,/</,{{q x }CI*I) m- l,{s x )}^& m- l,{s m }(i;D m- 1}) 

5 2 = Si{I 2 •-»• 100} = {/i i-> 10000, 7 2 ^ 100} 
/?2 = /?iU/mY(L 2 ,(0,S 2 ,/?i)) = {Ol 2 ^ (fi,ai)} 

L 3 = react ions (I 3 i-> 1,S 2 ) =0 

5 3 = S 2 {/ 3 ^1} = {/!^ 10000, 7 2 ^ 100, 7 3 ^ 1} 
R 3 =R 2 U init (L 3 ,(0,S 3 ,R 2 )) = R 2 

L4 = react ions (I4 1— >■ 1,S 3 ) =0 

5 4 = S 3 {/ 4 i-> 1} = {/1 i-> 10000,/ 2 ^ 100,/ 3 1-> l,/ 4 !-»• 1} 
/? 4 = /? 3 U/mY(L 4 ,(0,S 4 ,/? 3 )) =/? 3 
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with Ll = {{I 2 ' y l,h ^ IhWM,^}^ ^ h{®m z J h- >-> 1}) 

a i = /?ro pe/wity ( Ll , S 2 ) = 1 0000000 
fl =0 + rfe/ay(10, 10000000). 

Now, the reaction O in <\P\) X ^> T is This means that (jPD A ^> T is derived by means of an 
application of the ( Reaction rule[ ) as follows, where E = x,y,z,w, Si = {/ 2 i-» Ml ^ 1}> and S 2 = 
{{^}dp ^ l,{®4}a*5 ^ l,{s ffl }a^D m. 1}: 

((5i,10,frf,5 2 ),10,fi) =nexf(0,5 4 ,/g4) (£' h (0,S 5 , 7? 5 )) = h (0,S 4 ,/? 4 ),Si) 

£ h (0, S 4 , R4) 1Q ' (S ' ,10 '^' 52) > «orm(£' U /n(S 2 ) h (S 2 © ((f 1 , S s , R s ) G Si ) ) ) 

where 

S 5 = {h ^ l,{a,,.s m }<$l> 1 ^ 9999,/ 2 ^ s x }«JI> ^ 99,/ 3 ^ 1, 

7 4 ^ l,{® t }«i» ^ ^ 1} 

^5 = {<?l 2 i-)- (fi,ai),Oi 1 y (t 2 ,a 2 ),0 2 h-> () 3 ,a 3 ),6> 3 h-> (74,04)} 

with 

01 = ({{»„.o m }^,I)^l,/ 2 ^l},10,W,{{s x }(Ip^l,{s x }^,I)^l ! {s m }(I^I)^l}) 

02 = ({s x (s x ),s x }<IJI> -> M h- l},10,irf,{{s x }(IJD ^ l,{s x }< m D m- l,{o«}< D ^ 1}) 

03 = ^ l.l^M^najD ^ 1}, 10,W,{{» X }<$1> h- l,K;}<£l> -> l,{» m }aj>}) 
a 2 = propensity{Oi,S 5 ) = 99990 

a 3 = propensity(0 2 ,S$) = 990 
a 4 = propensity(0 3 ,S 5 ) = 989901 

f 2 = + </e/ffy( 10,99990) 

f 3 = + <ie/a.y(10,990) 

f 4 = + Je/aj(10,989901) 

We can now compute the multiset of the new machine state: 

s 2 ©((? 1 ,s 5 ,/? 5 )esi) = s 2 ©((fi,s 5 ,/? 5 )e{/ 2 ^i,/i^ 1}) 
= s 2 ©((/i,s 6 ,/? 6 })e{/i^i}) 
= ^ 1 ^ i,{s m }a?i> ^ i}©(?i,s 7 ,« 7 ) 

= (h,S&,R%) 



with 










= propensity (0 1, Se) = 


s 8 


= S 6 U{{^}0 M»m}«w» ^ M»m}«P -> 1} 




= fi + (a 2 /a 2 )(? 2 -?i) 




= {Oi^(44)} 


/ 


= propensity (0 2 ,Sj) = 




= R 6 U{0 2 ^(t' 3 ,a' 3 )} 




= fl + 3 M)(*3-fl) 


#8 


=R 7 u{{{v ± m WJ h. i,{»,„}a;rD ^ i},5,/,{8W ^ 1} 


s 6 


= S 5 \{/ 2 ^1} 






s 7 


= S 6 \{/i^l} 


/ 


= Ar.r[y:=z] 
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6 Conclusions 

In this paper, we have presented an abstract machine for the Stochastic Brane Calculus. Instead of defin- 
ing an ad hoc machine, we have adopted the generic abstract machine for stochastic calculi (GSAM) re- 
cently introduced by Lakin, Pauleve and Phillips. According to the encoding technique we have adopted, 
membranes are flattened into a set of species, where the hierarchical structure is represented by means 
of names. In order to keep track of these names, and for dealing efficiently with multiple copies of the 
same species, we have introduced a new generic abstract machine, called COWGSAM, which extends 
the GSAM with a name environment and a copy-on-write optimization strategy. We have proved that 
the implementation of the Stochastic Brane Calculus in COWGSAM is adequate with respect to the 
stochastic structural operational semantics of the calculus given in Q . 

We think that COWGSAM can be used for implementing other stochastic calculi dealing with nested 
structures, also beyond the models for systems biology. In particular, it is interesting to apply this ap- 
proach to Stochastic Bigraphs (5), a general meta-model well-suited for representing a range of stochas- 
tic systems with compartments; in this way we would obtain a General Stochastic Bi graphical Machine, 
which could be instantiated to any given stochastic bigraphic reactive system. However, such a machine 
would not scale well, as in general the COW strategy may be not very useful; thus, we can restrict our 
attentions to smaller subsets of BRSs, specifically designed to some application domain. For biological 
applications, the bigraphic reactive systems considered in [ 1 , 7 ] might be a more reasonable target. 

Another interesting question is about the expressive power of GSAM and COWGSAM. We think 
that GSAM correspond to stochastic (multiset) Petri nets, but COWGSAM could go further thanks to 
the possibility of creating unlimited new names during execution. Further work include comparison with 
other stochastic simulation tools dealing with compartments, like BioPEPA 151 . 
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A Some measure theory 

Given a set M, a family £ of subsets of M is called a o-algebra if it contains M and is closed under 
complements and (infinite) countable unions: 

1. Mel; 

2. A G £ implies A c G E, where A c =M\A; 

3. {A,} iGN C £ implies UieN^i £ s - 

Since M G £ and M c = 0, G £, hence £ is nonempty by definition. A a-algebra is closed under countable 
set-theoretic operations: is closed under finite unions (A,Be£ implies AU5 = AU5U0U0U--- G £), 
countable intersections (by DeMorgan's law AHB = (A c U B c ) c in its finite and inifite version), and 
countable subtractions (A, B G £ implies A \ B = A n B c ' G £). 

Definition A.l (Measurable Space). Given a set M and a o-algebra on M, the tuple (M, £) is called a 
measurable space, the elements ofL measurable sets, and M the support-set. 

A set Q. C 2 M is a generator for the a -algebra £ on M if £ is the closure of Q. under complement 
and countable union; we write a (£2) = £ and say that £ is generated by Q.. Note that the a-algebra gen- 
erated by a Q. is also the smallest a-algebra containing Q., that is, the intersection of all a-algebras that 
contain Q.. In particular it holds that a completely arbitrary intersection of a-algebras is a a-algebra. A 
a-algebra generated by £2, denoted by a (£2), is minimal in the sense that if £2 C £ and £ is a a-algebra, 
then a (£2) C £. If £2 is a a-algebra then obviously a (£2) = £2; if £2 is empty or £2 = {0}, or £2 = {M}, 
then a(£2) = {0,M}; if £2 C £ and £ is a a-algebra, then a(£2) C £. A generator £2 for £ is a base fori, 
if it has disjoin elements. Note that if £2 is a base for £, all measurable sets in £ can be decomposed into 
countable unions of elements in £2. 

A measure on a measurable space (M,£) is a function /!:£—)• K+, where R+ denotes the extended 
positive real line, such that 

1. JU(0)=O; 

2. for any disjoint sequence {A^} ;e / C £ with / C N, it holds 



M(Ute/tf/) = Ei€/M(tf/)- 



M. Miculan and I. Sambarino 



99 



The triple (M,E, ju) is called a measure space. A measure space (M,E,/i) is called finite if /x(M) is a 
finite real number; it is called o -finite if M can be decomposed into a countable union of measurable sets 
of finite measure, that is, M = U/er^i' f° r some / C N and /I (iV;) E M + for each i E /. A set in a measure 
space has o-finite measure if it is a countable union of sets with finite measure. Specifying a measure 
includes specifying its domain. If /I is a measure on a measurable space (M,E) and E' is a a-algebra 
contained in E, then the restriction // of /I to E' is also a measure, and in particular a measure on (M',E'), 
for some M'CM such that E' is a a-algebra on M'. 

Given two measurable spaces and measures on them, one can obtain the product measurable space 
and the product measure on that space. Let (Mi,Ei) and (M 2 ,E2) be measurable spaces, and pL\ and 
H2 be measures on these spaces. Denote by Ei (g> E2 the a-algebra on the cartesian product Mi x M 2 
generated by subsets of the form By x fi 2 , said rectangles, where B\ E Ei and B2 E E2. The product 
measure pL\ ® /i 2 is defined to be the unique measure on the measurable space (Mi x M 2 ,Ei ®Ea) such 
that, for all By E Ei and B 2 E E 2 

(/ii(8)/X2)(fli xB 2 ) = /Xi(fli)-ju 2 (52) 

The existence of this measure is guaranteed by the Hahn-Kolmogorov theorem. The uniqueness of the 
product measure is guaranteed only in the case that both (Mi,Ei,jUi) and (M2,E2,;U 2 ) are a-finite. 

Let A(M,E) be the family of measures on (M,E). It can be organized as a measurable space by 
considering the a-algebra generated by the sets {/I E A(M,E) : ju(5) > r}, for arbitrary 5 E E and r > 0. 

Given two measurable spaces (M,E) and (N,&) a mapping /: M — >• AT is measurable if for any 
T E 0, f~ l (T) E E. Measurable functions are closed under composition: given f:M—>N and g: N ^ O 
measurable functions then g o / : M — >• is also measurable. 



B Proof of Prop. 5.4 



Let P = Q^.rlTodQn.alaodP'D °P"D; then, (\P\) X = £ h species^(P) (0,0,0), where 

species 9tX (P) = {s(^.T\To)<lyt>}Uspecies {x}ty (^ n .a\o Q (lP'\) ° P") 
= {({^.t}U^(t ))(!J;D h-> l}UjpedeJ W!y (o„.a|a (I/ >/ D)Ujpedej W)) ,(P") 
= {({^.t}U^(t ))(!J;D ^ 1}U {*(»„. a|a )«>;i) ^ l}Uj/7edej {Xiy})W (P')Ujpeciej Wiy (P") 
= {({o^.t}Uj(to))CIJD ^ l}U{(j(s„.a)Uj(a ))^5 ^ l}U^ed^ 3 , }jW (P , )U^c^ {4iV (P ,/ ) 
= {({^.z}Us(z ))d x y D 1 ^ 1}U {({»„. a} U^(a ))(I 3 v ;D i-> l}Ujpede Jw})W (P')U^edej {4iy (P") 

Let/i = ({si.T}Uj(T ))ap,/ 2 = ({0 B .ff}Uj(0b))O,/j*=Jpe«M^ 

then dpo, = h ■-> 1 e/2 !->■ 1 e/p* e/p" e (o, 0,0) = / 2 ^ 1 ®j p ®Jp» e (0,5',/?') = j p > @j p » e (0, 5",/?") 

where 

L' = reactions(I\ i->- 1,0) = 

5' = 5{/i H- 1} 

/?' = init(L', (0,5', 0)) = 

L = reactions^ > 1,5') 

= ({/ 2 ^ l,/i ' ^ l},r n ,/,{(j(T)Uj(T )Uj(a)Uj(a ))(ip^l}) 
5" = 5'{7 2 ^1} = {/ 1 ^1,/ 2 ^1} 
/?" = inh(L",(0,S",R')) = {0 L (-»• (t l9 ai)} 
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with L = ({h i-> l,h ^ l},r n ,f,{(s(z)Us(z Q )Us(a)Us(a Q ))(l x y V i-> 1}). Now, the reaction O in 



d^Dx — T is C?l. This means that <\P\) X — T is derived by means of an application of the ^Reaction rulel) 



as follows, where S { = {h ^ l,h ^ 1}> $2 = {0(t)Uj(to) Uj(cj)U s(o ))4*i i-> 1} and f = T[w:=x] 

((Si,r n ,f,S 2 ),ai,ti)=next(0,S",R") (E 1 h (0,S m ,R m )) = cow(E h (0,S" X)^^ 
E h (0,S",P") ai ' (5l ' r "' / ' 52) > norfft (^u/n(5 2 ) h (/(5 2 e((fi,5 /,, ,/? ,,, )e5i)))) 
where S'" = S" and R'" = R". 

(f(s 2 e((t l ,s"',R"')es l )))) = 

= /((*(*) u^(to) Ui(a) u*(cto)^d) e ((?i,s"',/0 e {/j ^ i,/ 2 ^ i}) 

= /((*(T)U*(T )U*(ff)U*(ff )<ip)e ((?!,{/! ^0,/ 2 ^l},{O L ^(? 2 ,a2)})e{/2^1}) 

= /(0(t) u j(Tb) Ui(ff) u j(ob)<^i>) e (fi, {/i ^ o,/ 2 1 y 0},{O L i ^ (f 3 ,a 3 )}) 

= 0(t) U j(to) Uj(ff) U j(ob)^) © (fi, Ui ^ 0,/ 2 i ^ 0}, {<9 L i ^ (f 3 ,«3)}) 
= /©(f 1) {/iM.0,/ 2 M-0},{O i ^(/3,a3)}) 
Now let us define Q as Q = a\a \z\zo(iP'"i) °P', then (|0. v = species%^(Q) © (0,0,0) where 

species@j(Q) = speciesQ yX (o\oo\T\To<lP"D ° 

= specie ]X (o\oo\T\T (lP"D)Lispecies^(P') 

= {j(ct|(7o|t|to)^I)) h-> l}Ujpede JWiy (P")U^ecjej 0iX (P') 

= {s(T)Us(T )Us(a)Us(ob)<I*I> i-> \}UJp»UJ P > =J'UJp»UJp> 

And hence g = {J P i © / P « © J' © (/i , {/i i-> 0,/ 2 h-> 0},{Ol h-> (f3,a3)})]*. It remains to prove that r n = 
Hid([Q))- L et us notice that the derivation of P — > p. is actually as follows: 

, (loc) 

Sj.T->[sa T (comp) -> /X / ®g- g| ° b ([^]g B ..c<Dc b A*") P" -> ^ 

^ ^.T|T -> [»i] t8 ij,|l B„.CT|Ob<lP'l) °P" -> (At'@f"' <T|CTO ([»„]a B „. t 7®a Ai")) B „.a|cTo(!/"l)®/'"M'" 

(loc) 



s^.T|T (I»„.a|(Tot!P'D °P"D -> v 
where v = ((M'<4"- CT| ^[Q„] CTO ,,.a® ff0 M'0)^ 

Hi = (iU / @p/"' CT|CTC ([» ;1 ]c7S„. C 7®C7oJ U// ))s„. C 7|c7o<!P'D©P"i U/ " *nd jU 2 = [^]t^.T®T M- Then: 

v w ([a|a |T|T ap"D oP]) = ( Ml ©g^,^,, ^(MobM^'* "^D 
= Aii iy ([a|a |T|T (!P ,, Dop']) + Atlern ([a|ao] x [P] x [P']) -^([tIto])/^ 

= jU leXn ([(T|(To] X [P] X [P'])-Ai2 B(i ([T|To])A„ 

= (0* , e^ l< *([oJa O ..o<Do b M'')U([ff|ob] x [P'] x [P']) + <,([a|a ] x [P] x [P'])) -^([tKo])/^ 

= (([Oja B ,,a®a M"K(M<*>]) +<,([ff|Ob] X [P'] X [P'])) • ([^] TB|i . T ® To /i)^ ([t| T ] )/r„ 
= (([0„]ak([(X|C7o]) + <([<y^ 

= (r„ + <([cr|ob])+^(Wob] x [P'] x [P"]))-(r„ + M Si (M^o]))A, I = r n 

where the last equivalence holds because /4' n ([oiob]) = At^„([ CJ l a 'o] X [P'] x [P"]) = /i S( l([t|t ]) = 
because we assumed that the reaction does not involve either Ob nor To- 



